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ABSTRACT 

We discuss the effects of finite sky coverage and the uncertainty in extracting information 
about the power spectrum from experiments on small angular scales. In general the cosmic 
variance is enhanced by a factor of An /A, where A is the solid angle sampled by the exper- 
iment. As a rough guide, an experiment with sensitivity peaking at the ith multipole has 
to cover > i independent patches to have a lower "sample variance" than for a whole-sky 
measurement of the quadrupole. Our approach gives a relatively simple way of attaching 
an error bar to the theoretical prediction for a particular experiment, and thereby com- 
paring theories with experimental results, without the need for computationally-intensive 
Maximum Likelihood or Monte Carlo calculations. 

Subject headings: cosmology: cosmic background radiation — methods: statistical 
Introduction 

Since the discovery of large angular scale anisotropies in the microwave background (Smoot 
et al. 1992) effort has concentrated on constraining the primordial power spectrum of 
fluctuations using the results from degree-scale experiments (e.g. Gorski 1992, Vittorio & 
Silk 1992, Gouda & Sugiyama 1992, Dodelson & Jubas 1993, Gorski et al. 1993, White 
et al. 1993, Crittenden et al. 1993). It has become apparent (e.g. Gould 1993, White et 
al. 1993) that the "cosmic variance" , due to sampling only one universe, is a fundamental 
limitation on the utility of measurements on the largest scales. This uncertainty comes 
about because the fluctuations we observe are a single realization of a random variable. 
The theory can predict properties of the ensemble, but not the individual realizations. For 
example the quadrupole is drawn from a distribution with only 5 degrees of freedom, 
so that the purely theoretical error on the measured r.m.s. quadrupole is large. 

Usually it is assumed that the "cosmic variance" is negligible for small scale (< 1°) ex- 
periments which are sensitive only to high multipoles in the AT/T expansion. However, 
this assumption is only true for a whole-sky measurement. It is clear that the variance 
from only sampling a fraction of the sky will be larger. We feel that the importance of this 
"sample variance" has not been emphasized until now. 

In this letter we address the issue of the effects of finite sky coverage for small scale 
experiments and show that while in principle the sample variance can be small, currently 
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it is a major source of uncertainty. We find that the naive expectation is essentially correct, 
namely that the sample variance cr^g^^ is related to the cosmic variance a^^^ via 

^sam - (47rM)cT2os , (1) 

where A is the solid angle covered by the experiment. 

Consider an (idealized) experiment which is sensitive only to the £th. multipole. For a 
Harrison-Zel'dovich spectrum of fluctuations, in general we have (Tcos,i/ fJ-i ~ (•^+ V2)~^^^, 
where is the theoretical r.m.s. value of the £th multipole. If this experiment covers N 
independent patches, each of size ~ steradians, then we have 



'-''sam,^ 
fJ't 




Vm- (2) 



Recall that the cosmic variance for a (whole sky) measurement of the quadrupole gives 
an uncertainty of order the mean. Thus, obtaining a measurement with a theoretical 
uncertainty below that of the full-sky quadrupole requires N > i. 

Since this rule of thumb is not satisfied for most of the current small angular-scale experi- 
ments, we were motivated to investigate the effect of the "sample variance" in more detail. 
We find for the MAX (Meinhold et al. 1993, Gundersen et al 1993) and SP91/ACME 
(Gaier et al. 1992, Schuster et al. 1993) experiments a theoretical uncertainty of approx- 
imately 25-30% in the prediction of AT/T on the appropriate scale. Many authors use 
a Bayesian technique (e.g. Bond et al. 1991), or a Monte-Carlo simulation (e.g. Gorski et 
al. 1993), in analyzing these experiments. Such techniques automatically account for the 
sample variance, however sometimes the (dominant) source of the uncertainty is obscured. 
In the language of Bayesian analysis, our sampling uncertainty indicates that the data from 
small patches of the sky are highly correlated (the amount of experimental "information" 
is low), which leads to sensitivity to the "prior distribution" used in the analysis. 

We should stress that this sample variance has nothing to do with the experimental preci- 
sion of the measurements, but simply the fact that they may not necessarily have covered 
enough of the sky to provide a good estimate of the r.m.s. value of AT/T on the relevant 
angular scale. In the rest of this letter we develop this argument more rigorously, and give 
estimates for specific experimental configurations. 

Derivation of the sample variance 

Consider the correlation function at zero-lag, or [(AT/T)„ns]^, as estimated by measure- 
ments over a patch of the sky subtending a solid angle A and which we shall call A^. We 
include the effects of finite beam size and possible "chopping" of the beam in our definition. 

The theory predicts a probability distribution for A^, but not its value in our universe. 
Write Co as the mean correlation function at zero-lag: (A^) = Cq, where the angled 
brackets represent an average over an ensemble of "universes". This is a measure of the 
amplitude of the theoretical power spectrum at the scales probed by the experiment. 
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Under the assumption that the original fluctuations are gaussian, it is straightforward to 
calculate 

" 21 (3) 

■T(ni)T(ni) ) (T(n2)T(n2) ) + 2 (T(ni)T(n2) ' 



^ / dnidn2 

I A 



Here T{n) refers to the temperature difi'erence assigned to the direction h by the exper- 
iment, including finite beam width and possible beam "chopping". The angled brackets 
again represent an average over an ensemble of "universes" . To fix the size of the uncer- 
tainty in translating from the experimentally measured to Co we calculate the sample 
variance, which is simply the second moment about the mean for A^: 

= a(J) = (Ai> - (A^)^ 



[ dVLidVL2C'^{fii-fi2) 

J A 

where C(cos6') is the 2-point correlation function 

C(ni • n2) = (T{nr)f{h2) 



(5) 



47r „ 
e=2 



Here a| = \ciem\ , where the aim ai'e the multipole coefficients of the temperature 
fluctuation on the sky and is the window function appropriate to the experiment under 
consideration (see Bond et al. 1992, Dodelson Sz Jubas 1993, White et al. 1993). 

One can also derive the third moment about the mean: 

A^^^ = 13 / d^id^2dn3 C{fii ■ fi2) C{h2 ■ hs) C{fi3 ■ ni) . (6) 

In general the distribution for Aa will be positively skewed. Similar expressions can be 
written for all the higher moments. 

Simple Models 

For the purposes of illustration we consider the case of a gaussian autocorrelation function 
(GACF). Specifically we take 

C(cos6') = Coexp[(cos6' - l)/el] ~ Co exp [-6* 726*2] , (7) 

and we assume 9c <^ l- The correlation function, C(cos^), for both the SP91/ACME 
and MAX experiments (to choose specific examples) with CDM-like power spectra can be 



adequately approximated by equation (7). We should emphasize that our C(cos6') has 
already been convolved with the telescope beam, thus 9c is a (calculable) function of the 
power spectrum and the experimental parameters. Often the unconvolved C {cos 9) is taken 
to be a GACF, which is not a good assumption. 

We see from equation (4) that we expect A^-* ~ Cq9c/A for 9^- To explore this further 
consider a circular region defined by < ^ < 6*0 and < < 27r, so that A = 27r(l — cos ^o)- 
Using equations (4) and (7) we find that 

= ^^ ^^\\2 f dc,dc2 exp[2(ciC2 - l)/9l]lo{2s^S2/9l) , (8) 

(i -cos 6^0 j JcosSo 

where q and Si are cos 9i and sin 9i respectively, and Iq is the modified Bessel function of 
order zero. We have the two limiting cases 



^^I^ - i 1 /4.\ . . (9) 




We show A^"* and \{A.'k / A)Cq9'1 for 9c = 1° in figure 1. The expected result is recovered 
in the limit of one measurement: for a gaussian random variable x, the variance of is 
twice the square of the mean of x'^. Notice that above a few degrees the second limiting 
form is a good approximation to the exact result. This means that after a few correlation 
lengths the sample variance, A^\ approaches the cosmic variance, ^Cq6'^, as An/A, so 
that there is always something to be gained by sampling more sky. 

While the region A defined above clearly illustrates the simple scaling of the sample vari- 
ance with sky coverage, it is not a good approximation to the scan strategy of SP91/ACME 
and MAX. To consider these experiments more properly we here make a simplified model 
of the scan strategies. In both cases we can approximate the scan pattern as a single line 
of (angular) length a. For our GACF assumption we find 



Ai^' = 2 r**CVos 1^.-^,1) 

Jo a a 

^2C„^^{v^.rf(£)-|[l-exp(-„V''2)]} 



(10) 



For small a we again recover the limit A*^^) = 2Cq while for a larger than a few 9^ the second 
term is small and the error function is approximately unity. Comparing with equation (9) 
we see that the effective solid angle A for this linear scan is v^^c x Putting in a MAX 



correlation angle 9c — V2° ^-nd a scan length of a = 6° we find that Osam = y A^"* ~ O.SCq; 

that is, the uncertainty is approximately 50% of the mean [(AT/T)i-ms] = {^a) = Co- In 
the case of SP91/ACME, 9c — 2° and a ~ 9 x 2°1, giving a relative error of approximately 
60%. The actual value of 9c depends on the precise cosmological model assumed, but the 
values used above are representative for CDM and these experiments. 
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The implication of this simple analysis is that the theoretical predictions for AT/T for 
the experiments considered should be assigned a relative error of approximately 25-30% 
due to the "sample variance". (However recall that the distribution of will be non- 
gaussian; see e.g. Scaramella & Vittorio 1991, Cayon et al. 1991, White et al. 1993.) 
Other experiments have similarly non-negligible sample variances because of incomplete 
sky coverage. This "range" of theoretical predictions should be considered in comparing 
theory with small scale experiments. Only COBE and MIT approach the unavoidable 
cosmic variance limit for the scales which they probe. Note, however, that the galactic cut 
made by the COBE team > 20°) increases the relative error on AT/T by 23% over 
what it would be for a whole sky measurement. 

Conclusions 

Our analysis makes clear that the potential gain in theoretical precision on smaller angular 
scales is in fact only realized once a significant fraction of the sky has been covered. Thus 
the effects of finite sky coverage are serious for all of the current small scale experiments: 
VLA, ATCA, OVRO, MSAM, MAX, SP91/ACME and (on larger scales) Tenerife. Includ- 
ing the "sample variance" in the theoretical predictions for these experiments may help 
to reconcile the apparently contradictory results being reported on small angular scales 
(Gaier et al. 1992, Schuster et al. 1993, Gundersen et al. 1993, Meinhold et al. 1993). Cur- 
rent analyses of the data (including Monte-Carlo and Bayesian analyses) correctly take 
into account the "sample variance" , but provide little physical insight or guide to the size 
of the effects. We hope that our approach is sufficiently straightforward that it will make 
the comparison of theory and measurement more transparent for these small scales, and 
perhaps also act as an aid for the design of future experiments. 
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Figure Caption 

(2) 

Figure: The sample variance, , from equation (8) as a function of sky coverage, 6*0, 
for a GACF with 0^ = 1° (solid line). Also shown is the naive expectation (47r/A)cr^os 
(dashed line) where A = 27r(l — cos^o) is the solid angle covered by the region sampled 
and (T^Qg = ^Cq9'^ is the cosmic variance. This is a good approximation to the full result 
above a few correlation lengths. 
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